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ABSTRACT 

Single degenerate model is the widely accepted progenitor model of Type la super- 
novae (SNe la), where a CO WD accretes hydrogen-rich material from its companion 
to increase its mass. The companion may be a main-sequence star or a subgiant star 
(WD -I- MS). When the CO WD approaches the Chandrasekhar mass limit, it ex- 
plodes as a SN la and a part of supernova ejecta collides into the companion envelope. 
After the impact of the ejecta, the companion survives and may show some special 
properties. A good way to verify the single degenerate model is to study the inter- 
action between supernova ejecta and its companion, and/or search companion in the 
rem nant of a SN la. 

iMeng. Chen fc Hani (|2009t ) comprehensively and systemically studied the WD -f 
MS system by detailed binary evolution calculations. Following their studies, we have 
carried out a series of binary population synthesis studies about the properties of the 
companions of SNe la for different metallicities Z. We present the distributions of the 
masses, Aff^, the radii, -Rf^, of companions, and the periods, Psn, and the ratios 
of separations to radii, A/flf^, of WD -I- MS systems for various Z at the moment 
of supernova explosion. These parameters can be applied to constrain the numerical 
simulation of the interaction between the ejecta of a supernova and its companion. 
We also show the distributions of some integral properties of companions, i.e. the 
mass, the space velocity and the surface gravity, for various Z after the interaction. 
The distributions may help to search companion in a supernova remnant. All the 
parameters above significantly change with Z . 

Incorporat ing the simulati o n resu lts o f interaction betwee n supernovae ejecta and 
companions in [Marietta et al.l ()2000D and iKasen et al.l (|2004| ) into our binary popu- 
lation synthesis study, we found that more than 75% of all supernovae have enough 
polarization signal which can be detected by spectropolarimetric observations. We 
also found that 13 to 14 per cent SNe la belong to the supernovae like 1991T, which 
is consistent with observations within errors. This may indicate that SNe 1991T-like 
have not any special properties in physics except for the viewing angle of an observer. 

Key words: binaries: close-stars: evolution-supernovae: general-white dwarf: metal- 
licity 



1 INTRODUCTION 

Type la supernovae (SNe la) play an important role in 
astrophysics, especially in cosmology. They appear to be 
good cosmological distance indicators and are successfully 
ap plied to determine cosmological parame ters (e.g. Q and 
A: lRiess et aLlll998l : IPerlmutter et al.lll999l ). There is a lin- 
ear relation between the absolute magnitude of SNe la and 
the magnitude difference from maximum to 15 days after B 



maximum ligh t. The relation is known as Phillips relation 
(iPhillipj|l99a ) and adopted when SNe la were taken as the 
distance indicators. In this case, Phillips relation is assumed 
to be valid at high redshift, although it was obtained from a 
low-redshift sample. This assumption is precarious since the 
ex act nature about SNe la is sti l l unclear (see the reviews 
by iHillebrandt fc Nieme"ve3l2000l : lLeibundgutjl2000D . If the 
properties of SNe la evolve with redshift, the results for cos- 
mology might be different. Since metallicity may represents 
redshift to some extent, it is a good method to study the 
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properties of SN la at various redshift by finding tfie corre- 
lation between tlie properties and metallicity. 

It is widely believed that SN la is from the thermonu- 
clear runaway of a carbon-oxygen white dwarf (CO WD) 
in a binary system. The CO WD accretes material from its 
companion to increase its mass. When its mass reaches its 
maximum stable mass, it explodes as a thermonuclear run- 
away and almost h alf of the WD mass is converted into ra- 
dioactive nickel-56 (|Branchl 12004 1 . The mass of nickel-56 de- 
termines the maximum luminosity of SN la. The higher the 
mass of nicke l-56 is, the higher the maximum luminosity is 
l|Arnettll 19821 ). Some numerical and synthetical results have 
shown that metallcity may affect the fi nal amount of nickel- 
56, and thus the max i mum luminositv |Timmes et al.l l2003l : 
iTravaglio et~aLl l2005l : iPodsiadlowski et al.l |2006| ). There is 
also much evidence about the correlation between the prop- 
erties of SNe la and metallicity in o bservations (we list here 
some of the p apers on this stud y iBranch fc Berghl Il993l: 
Hamuv et all Il996|: IWang et al.1 1 19971 : ICappellaro et al] 



19971 : IShanks et al.ll2002i r 

Among all the suggested progenitor models of 
SNe la, the single-degenerate Chandrasekhar mode l 
l|Whelan fc Ibenll 19731 : iNomoto. Thielemann fc Yokoilll984l ) 
is the most widely accepted progenitor model at present. 
In the model, a CO WD accretes hydrogen-rich ma- 
terial from its companions until its mass reaches 
a mass ~ l.SYSMp, (close to Chandrasekhar mass, 
iNomoto. Thielemann fc Yokoil 1 1984 ). and then explodes 
as a SN la. The companion may be a main sequence 
star or a s ubgiant star (WD-I-MS ) or a red-giant star 
(WD -f -RG) llYungelson et al.1 Il995l: iLi fc vaii den Heir^ 
19971: iHachisu et al.l Il999allbl: INomoto et al.1 |l999l. l2003l: 



Langer et al. 120001: iHan fc Podsiadlowskil 12004 IChen fc Li 



2007*. '2009': Han' '2008* 'Me ng. Chen fc Hani |2009| : iLii et al 
2009; Wang ct al. 2009a, b). In this paper, we only focus 
on the WD -I- MS channel since it is the most widely 
accepted channel for SNe la (|Han fc Podsiadlowskil |2004| : 
iMeng. Chen fc Hani |2009| ). Much observational evidence 
shows the importance of the channel. For example, some 
WD -|- MS systems are sug g ested as the progenitor of SNe la 
iParthasarathv et al.|[200^ ). iHachisu fc Katol l|2003al lbl) sug- 
gested that supersoft X-ray sources (SSSs) may be good 
candidates for the progenitors of SNe la, and some of 
SSSs belong to WD-I-MS system. The discovery of the po- 
tential companion of Tycho's supernova seemed to verify 
the rehability of the WD -I- MS model (named Tycho G 
by Ruiz-Lapuente et al.l |2004 ) . Recently, iHernandez et all 
(|2009l ) stressed further the companion nature of Tycho G 
by analysing the chemical abundances of Tycho G. However, 
iFuhrmannI (|2005l ') argued that the first discovered "compan- 
ion" might be a Milky Way thick-disk star which is coinci- 
dentally passing the vicinit y of the remnant of Tycho's su- 
pernova, (ihara et al.1 l|2007h also argued that Tycho G may 
not be the companion of Tycho's supernova since the star 
did not show any special properties in its spectrum, which 
should be contami nated by supernova e.j e cta and show some 
special characters (iMarietta et al.ll2000l : iBranchI [20o3 ) . So, 
more evidence to confirmed the companion nature of Tycho 
G is needed. 

The knowledge about the companions of SNe la after 
explosions is still unclear. Generally, the supernova ejecta 
in the single degenerate model collides into the envelope of 



its companion and strips some hydrogen-rich material from 
the surface of the companion. After the collision, the com- 
panion gains a kick velocity, which is much smaller than 
orbital velocity, and leaves ex plosion center at a velocity 



orbital velocity, and leaves ex plosion center at a velocity 
similar to its orbital velocitv jcheng 1974: Wheeler et al 
' FrvxeU fc Arnettiri98ll:lTaam fc Frvxell.lll9 84':'Chugai 
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198e 



Liven. Tuchman fc Whelledl 19921 : [Langer e t al. 2000) 
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200Cf) ran several high-resolution two- 
dimensional numerical simulations of the collision between 
the ejecta and the companion, where the companion is a 
MS star, a subgiant (SG) star or a red giant (RG) star. 
They found that about 0.15 Mq - 0.17 Mq of hydrogen- 
rich materials are stripped from the surface of a MS or a 
SG companion and n a sense of the collision, there is no 
difference whatever the companion is a MS star or a SG 
star. The amount of the stripped material from red-giant 
companions is even more than O.SMp, i.e. mo r e than 96% 
of their envelopes is lost. iMeng. Chen fc Hani (|2007f ) used 
a simple analytic method but a more physical companion 
model to simulate the interaction and found that the min- 
imum value of the stripped material from a MS compan- 
ion is diminished from 0.15 Mq to 0.035 Mq. They sug- 
gested that the structure of the companion at the explo- 
sion moment may be an important factor to determine the 
mass of t he stripped material . The structure of the com- 
panion in iMeng. Chen fc Hani (|2007l) is obviously d i fferent 
from that of a solar model used in IMarietta et al.l (|2000l ) 
since the companion still do not reach thermal equilibrium 
because of mass transfer. The reduction of the stripped 
m aterial inlMeng^ Chen fc Hani (|2007l ) compared with that 
in IMarietta et al. (|2000t ) primarilv re sults from the pre- 
explos ion mass loss. The MS model in iMeng. Chen fc Hani 
(|2007h is from a more massive star experiencing a mass-loss 
phase before supernova explosion, which leads to a more 
compact companion star whose material is more difficult to 
strip t han it i s in a solar model a s used in IMarietta et al.l 
(|2000l ). Since iMeng. Chen fc Haiil (|2007l ) did not consider 
the thermal energy imparted by the ejecta into the com- 
panion envelope, which likely heats and vaporizes a part of 
the envelope and thereby increases the amount of stripped 
material, 0.035 Mq should be a conservative lower limit. If 
enough hydrogen-rich materials are stripped, they should re- 
veal themselves by narrow H a emission or absorption line i n 
later-time spectra of SNe la (IChugailll98^ : lFilippenkdll997l '). 
However, the Ha line was not detected and the amount of 
the stripped material was constrained to b e less than 0.0 2 
Mq by observations (|Mattila et al.1 l2005l : iLeonardI [20071 ). 
which is much smaller than 0.15 Mq - 0.17 and is even 
smaller than the lower limit obtained by'Meng, Chen fc Hanj 
( 2007.). In addition, the simula tion in iMarietta et al.. (2000!) 
and lMeng. Chen fc Hani l|2007l ) showed that the luminosity 
of supernova companion after collision by supernova ejecta 
should dramatically rise to a level much higher than that of 
Tycho G by about three orders of magnitude. 

Owing to the discussions above, a detailed numerical 
simulation on the interaction between supernova ejecta and 
its companion should be very important, while the subject 
is closely related with the prope rties of the secondar i es be - 
fore SNe la explosion. Recently, iMeng. Chen fc Hani (|2009[ ) 
(hereinafter Paper I) calculated a dense grids of binary sys- 
tems by detailed binary evolution for different metallicities 
and gave the initial parameters of the systems leading to SNe 



Companion Stars of Type la supernovae 3 



la in an orbital period-secondary mass (log Pi, A/2) plane. 
The aim of this paper is to, using the results in Paper I, 
show the properties of secondaries before and after the SNe 
la explosion, which can provide help to do detailed numerical 
simulations of the collision between supernova ejecta and its 
companion, and/or search the companion in the explosion 
remnant of a SN la. 

The paper is organized as follow. We simply show our 
binary population synthesis (BPS) method in section [2] and 
the BPS results in section |3l In section |4l we briefly discuss 
our results and summarize the main conclusions in section 
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Figure 1. The distribution of the final masses and the orbital 
periods at the moment of explosions. Common envelope ejec- 
tion efficiency acE = 1-0. The posi tion of a recurrent nova, U 
Sco, is ind icated by the filled square l lSchaefer fc Ringwald|[l995l : 
iHachisu et al.ii2000ai ibi). Top: Z = 0.03; Bottom: Z = 0.02. 
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Figure 2. Similar to Fig[T]but for Z - 

(Bottom). 



: 0.01 (Top) and Z = 0.004 



2 BINARY POPULATION SYNTHESIS 



2.1 



The parameters of WD+MS systems leading 
to SNe la 



Incorporating the prescription in lHachisu et al.l (Il999al ) on 
the accretion of the hydrogen-rich material from its com- 
panion onto a WD in t o th e Eggl cton's stellar evo l ution code 
(Ega ieton|[l97ll . Il97l [l973l ) , iMeng. Chen fc Hani (|2009l ) cal- 
culated more than 25,000 binary evolutions of WD-I-MS 
channel. In the channel, a companion fills its Roche lobe 
in MS stage or in Hertzprung gap (HG) and transfers some 
of its mass onto the surface of the WD. As a consequence, 
the mass of the WD incre ases gradually. If the mass of 
the W D reaches ~ 1.378Mq (|Nomoto. Thielemann fc Yokoil 
1984). they assume d that the WD explodes as a SN la. 
Me ng. Chen fc Hani (|2009l ) provided the initial parameter 
spaces of the WD -|- MS systems leading to SNe la in an or- 
bital period - secondary mass (logPi,M2) plane. Their cal- 
culations are comprehensive and systematical, and various 
properties of the companion stars with different metallicities 
were obtained but not sorted for publishing. In this paper. 
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Figure 3. Similar to Fig [T] but for Z 
0.0003 (Bottom). 



0.001 (Top) and Z 



Figure 4. Similar to Fig[T]but for Z = 0.0001. Top: acE = 1-0; 
Bottom: acE = 3.0. 



we extract the properties from the data files of the calcula- 
tions and incorpor ate them into the rapid bi nary evolution 
code developed bv lHurlev et all (|2000l . 1200 j ) to obtain the 
various properties of the companions at the moment of ex- 
plosion. 

2.2 Common envelope 

Common envelope (CE) is very important for the forma- 
tion of WD -I- MS systems. We firstly introduce the treat- 
ment for CE in this paper. Hereinafter, we use primordial 
to represent the stage before the formation of WD-I-MS sys- 
tems and initial for WD-I-MS systems. During binary evo- 
lution, the primordial mass ratio (primary to secondary) 
is crucial for the first mass transfer. If it is larger than 
a critical mass ratio, qc, the first mass transfe r is d ynam- 
ically unstable and a CE forms ( jPaczynnskij Il976l ). The 
ratio gc varies with the evolutionar y state of the primor- 



1987 




2002 





Webbinklll988l: lHa.n et"aLl 1200 j T IPodsiadlowski et all 
Chen fc Han " l2008l ). In this study, we adopt qc = 



4.0 when the primary is on MS or HG. T his value is sup- 
ported by detailed binary evolution studies (|Han et al ] |200G| : 
IChen fc Hani |2002| . |2003| ). If the primordial primary is on 
FGB or AGB, we use 



[1.67 -X 4-21 



1/2.13, 



(1) 



where is the core mass of primordial primary, and x = 
dlni?r/dlnMf is the mass-radius exponent of primordial 
primary and varies with composition. If the mass donors 
(primaries) are naked helium giant s, Qc = 0.748 based on 
equation (|T]) (see lHurlev et al ] |2002l for details). 

Embedded in the CE are the dense core of the primor- 
dial primary and the primordial secondary. Due to frictional 
drag with the envelope, the orbit of the embedded binary de- 
cays and a large part of the orbital energy released in the 
spira l-in process is injected into the envelope l|Livio fc Sokeij 
Here, we assume that the CE is ejected if 



(2) 



QCEAiJorb ^ l-Bbind 
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where ABorb is the orbital energy released, iSbind is the bind- 
ing energy of common envelope, and acE is CE ejection effi- 
ciency, i.e. the fraction of the released orbital energy used to 
eject the CE. Since the thermal energy in the envelope is not 
incorporate d into the binding energy, oce may be greater 
than 1 fsee lHan. Podsiadlowski fc Eggletonll 19951 for details 
about the thermal energy). In this paper, we set qce to 1.0 
or 3.0. 



2.3 Evolution channels 

There are three channels to produce WD -|- MS systems 
according to the situation of the primary in a primordial 
system at the onset of the first RLOF. 

Case 1 (He star channel): the primordial primary is in 
HG or on RGB at t he onset of the first RLOF (i. e. case B 
evolution defined by iKippenhahn fc WeigertI [l963) . In this 
CE is formed because of a large mass ratio or a 
convective envelope of the mass donor. After the CE ejection 
(if it occurs), the mass donor becomes a helium star and 
continues to evolve. The helium star likely fills its Roche 
lobe again after the central helium is exhausted. Since the 
mass donor is much less massive than before, this RLOF is 
dyn amically stable, resu l ting in a close CO WD-I-MS system 
fsee lNomoto et al.lil999l . |2003| for details). 

Case 2 (EAGB channel): the primordial primary is in 
early asymptotic giant branch stage (EAGB) (i.e. helium 
is exhausted in the core, while thermal pulses have not yet 
started). A CE is formed because of dynamically unstable 
mass transfer. After the CE is ejected, the orbit decays and 
the primordial primary becomes a helium red giant (HeRG) . 
The HeRG may fill its Roche lobe and start the second 
RLOF. Similar to the He star channel, this RLOF is sta- 
ble and produces WD -I- MS systems after RLOF. 

Case 3 (TPAGB channel): the primordial primary fills 
its Roche lobe at the thermal pulsing AGB (TPAGB) stage. 
Similar to the above two channels, a CE is formed during 
the RLOF. A CO WD + MS binary is produced after CE 
ejection. 

The WD -I- MS systems continue to evolve and the sec- 
ondaries may also fill their Roche lobes at a stage and Roche 
lobe overfiow (RLOF) starts. We assume that if the initial 
orbital period, -Porbi S'Hd the initial secondary mass, M2, of 
a WD + MS system locate in the appropriate regions in the 
(log P\M'2) plane for SNe la at the onset of RLOF (see Pa- 
per I), a SN la is then produced. The properties of the binary 
system at the moment of SN explosion are obtained by in- 
terpolation in the three-dimensional grid (M^di M2,\og P') 
of the more than 25,000 close WD binary system calculated 
in Paper I. 



2.4 Basic parameters in Monte Carlo simulation 

To investigate the statistical properties of companions before 
and after explosion, we carry out a series of binary popula- 
tion synthesis (BPS) studies for various Z by Hurley's rapid 
binary evolution code (Hurlcv ct al. 2000, 2002) incorporat- 
ing the results in Paper I. Since the code is only valid for 
Z ^ 0.03, only seven metallicities, i.e. Z = 0.03, 0.02, 0.01, 
0.004, 0.001, 0.0003 and 0.0001, are examined here. In each 
BPS study, lO'^ binaries are generated by Monte Carlo sim- 



ulation and a circular orbit is assumed for all binaries. The 

basic parameters for the simulations are as f ollows. 

(i) The initial mass function (IFM) of iMiUer fc Scald 

(|l979l ) is adopted. The p rimordial prima r y is g enerated ac- 
cording to the formula of lEggleton et al 



ir y IS g ei 



Mf = 



0.19X 



(1 _X)0-75 + 0.032(1 -X)0'- 



(3) 



where X is a random number in the range [0,1] and Mf is 
the mass of the primordial primary, which ranges from 0.1 
Mq to 100 Mq. 

(ii) The mass ratio of the primordial components, q, 
is a very important parameter for binary evolution while 
its distribution is quite controversial. For simplicity, we 
take a uniform mass-rat io distribution l|Mazeh et al.l 19921 : 
iGoldberg fc Mazeh|[l99i ): 



n(g) = 1, 



< (/ < 1, 



(4) 



where q = M^/Mf. 

(iii) We assume that all stars are members of binary 
systems and that the distribution of separations is constant 
in log a for wide binaries and falls off smoothly at close sep- 
aration: 



an{a) 



ao < a < ai, 



(5) 



where Qscp ~ 0.070, ao = IORq, ai = 5.75 x 10® 7?© = 0.13pc 
and m « 1.2. This distribution implies that the num- 
bers of wide binary system per logarithmic interval are 
equal, and that approximately 50% of the stellar systems 
are binary systems with orbital perio ds less than 100 yr 
(|Han. Podsiadlowski fc Eggleton|[l995l ). 

(iv)We simply assume a constant star formation rate 
(SFR) over last 15 Gyr. 



3 BINARY POPULATION SYNTHESIS 
RESULTS 

3.1 the companion masses and the orbital periods 
of WD + MS systems at the moment of 
explosions 

The secondary masses and the periods of WD -I- MS systems 
at the moment of explosion are basic input parameters when 
one simulates the interaction between supernova ejecta and 
its companion. We show the distributions of the masses and 
the periods for different metallicities Z in Figs.[T]to|31 where 
MwD = 1.378AfQ. Here, we only show the cases of acE = 
1.0 except for Z = 0.0001 (see the bottom panel in Fig. 
21 . For other Z, qce does not significantly affect the final 
distributions. 

From these figures, we see that generally, the compan- 
ions with a high metallicity Z have larger final masses, and 
roughly when Z > 0.004, the final periods center on 1 to 
1.6 days, smaller than those of the cases of Z ^ 0.004 (cen- 
tering on 2.5 to 6.5 days). These phenomena are mainly 
derived from the distributions of the initial parameters of 
WD -l-MS systems with various Z. Paper I has shown that 
a high metallicity leads to a more massive initial secondary 
mass(see Figs. 4 and 10 in Paper I), and the systems with 
Z > 0.004 have a shorter initial period than those with 
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Z ^ 0.004 (see Fig. 11 in Paper I). Meanwhile, metallicity 
also affects mass transfer pro cess. A high metalli city results 
in a lower mass transfer rate ijLanger et al.l[2000 ^. and then 
less material loses fro m binary system as optically thick wind 
l|Hachisu et al.lll996l ). which means that more transferred 
materials are accumulated on the WD. Therefore, mass-ratio 
inversion for a WD + MS system with a high metallicity is 
delayed, and then the WD + MS system is more likely to de- 
crease its orbital period, whi ch may lead to a relatively low 
period (see also Paper I and lMeng. Yang fc Genj|2009l ). 

The bottom panel in Fig. |4] shows the case for Z — 
0.0001 and acE = 3.0. In the figure, there are two "islands" 
in the (Mf^,logPsN) plane, i.e. a low-mass low-period one 
and a high-mass high-period one, which is significantly dif- 
ferent from the top panel in Fig. U We have not found the 
two islands in other cases. As shown in subsection 12.21 and 
paper I, there are three channels leading to close WD -I- MS 
binary systems, i.e. He star channel, EAGB channel and 
TPAGB channel. Metallicity and acE may systemically af- 
fect the channels (see Paper I for the influence of Metallicity 
and acE on the channels for details). A low metallicity and 
a high QcE, i.e. Z = 0.0001 and acE ~ 3.0, may result in 
double peaks in the distributions of the initial masses of CO 
WDs and the initial periods of WD -I- MS systems leading 
to SNe la (see Figs. 9 and 11 in Paper I). The low-mass 
peak and the low-period peak are from the EAGB channel, 
while the high-mass peak and the high-period peak are from 
the TPAGB channel. It is obvious that the low-mass peak 
of CO WDs and the low-period peak may lead to the low- 
mass low-period island in the bottom panel in Fig. |4] since 
more materials are needed to transfer onto CO WD from 
mass donor if a low-mass WD is adopted. For a similar rea- 
son, the high-mass high-period island in the bottom panel 
in Fig. m results from the high-mass peak and high-period 
peak from TPAGB channel. 
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When Z > 0.0001 and qce = 3.0, no WD + MS sys- 
tems are from TPAGB chann el, which results in the disap- 
pearance of double peaks (see lMeng. Chen fc H an 2009 for 
details). This can naturally explain why the double islands 
in Fig. [1] to the top panel in Fig. |4] disappear. 

These figures may also help to check whether some 
binary system observed can explode as SNe la or not. In 
these fi gures, a recurrent nova, U Sco, is indicated by fil led 
square (ISchaefer fc Ringwald|[l995l : iHachisu et a"l]|2000al lbh. 
The WD mass of U Sco is about 1.3 7Mq, and its compan ion 
is a MS sequence star of 1.5Mq (iHachisu et all l2000allbh. 
The o rbital period i s 1.230 56 days ( Schaefer fc Ringwaldl 
119951 ). IHachisu et all l|2000al lbl) studied the system carefully 
and concluded that the WD mass can grow until an SN la 
explosion is triggered after ~ 10^ yr. If the WD of U Sco 
may explode as a SN la finally, the companion should have 
a slightly smaller mass. Its period will decrease firstly, and 
then may increase if mass-ratio reverses. Then, its final po- 
sition in (Mf^, log Psn) plane will move to a lower mass 
from its present position at least, and may enters into the 
most probable area (see the bottom panel in Fig[T]). From 
its present position in (Mf^, log Psn) p lane, it is very likely 
for U Sco to explode as a S N la (see also IHachisu et al" ] |20oA 
iMeng. Yang fc Gendiioogl ). 



Figure 5. The distribution of the radii and the ratios of sepa- 
rations to radii of companions at the moment of supernova ex- 
plosion for the case of acE = 1.0. Cross represents Tycho G, 
which is a potential candidate of the companion of Tycho's su- 
pernova l lRuiz-Lapuente et al.ll2004l : lBranchll2004) . The length of 
the cross represents observational error. Solar symbol and star 
represent the main-sequence and subgiant companion model, re- 
spectively, which are applied to simul ate the collision betwe en su- 
pernova ejecta and its companion in [Marietta et ah I ll200[t) . The 
positio n of a recurrent nova, U Sco, is indicated by t he filled 
square fechaef er fc Rinewaldlll995lMHachisu et al.ll2000al lbl). Top: 
Z = 0.03; Bottom: Z = 0.02. 



3.2 Radii and the ratios of separations to radii of 
companions at the moment of explosions 

The radius of a companion, R^^ , and the ratio of sep- 
aration to the companion's radius, A/R2^, at the mo- 
ment of explosion are important parameters for simulating 
the inter action betw een supernova ejecta and the compan- 
ion (Mari etta et al., 2000). There is even a linear depen- 
dence of kick ve locity and the mass of stripped material 
on log(^/7?f ^) ([Marietta et all l2000l : iMeng. Chen fc Hani 
l2007f ). After the interaction, the companion reestablishes dy- 
namical equilibrium quickly while it is still in a process into 
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Figure 6. Similar to Fig[5]but for Z = 0.01 (Top) and Z = 0.004 
(Bottom). 



Figure 7. Similar to Fig |5] but for Z = 0.001 (Top) and Z 
0.0003 (Bottom). 



thermal equilibrium. The process into thermal equilibrium 
may last for 10^ - 10^ yr. It is similar that a companion 
star is also not back into thermal equilibrium at the mo- 
ment of supernova explosion since mass transfer is process- 
ing before supernova explosion. So, considering the amount 
of the stripped material from companions is small (see sec- 
tionUor next subsection, iMattila et al.|[2005l : iLeonar dl l2007l : 
iMeng^ Chen fc HanllioOTl ). the radius of a companion at the 
moment of supernova explosion may represent its radius af- 
ter the impact of a SN la to some extent, and then our 
results can also compare with observation directly. Under 
the assumption that the r adius of a c ompanion is equal to 
its Roche lobe radius (lEggletonl Il983l ) . we show the distri- 
butions of i?f^ and A/Rf^ at the moment of explosion for 
various Z and acE = 1-0 in Figs. [5] to |8] The cases with 
acB = 3.0 are similar to that with acE ~ 1-0 except for 
Z = 0.0001 (the bottom panel in in Figl}. 

The age of Tycho's supenrova is about 440 yr, which 
means that the suggested companion star of Tycho's supen- 
rova, Tycho G, is still not back into thermal equilibrium 
and its radius at present might be similar to that at the mo- 



ment of supernova explosion, also similar to its Roche lobe 
radius. Assuming that the radius of Tycho G at present is 
equal to its Roche lobe radius at the m oment o f super nova 
explosion, and using the equation in lEggletonl (|l983l ). we 
may obtain the of Tycho G. The cross in the figures 

represents the Tycho G (|Ruiz-Lapuente et al.|[20o3 : iBranchI 
[20043). The solar symbol and open pentacle denote the MS 
(a solar model) and SG companion model, which are ap- 
plied by [Marietta et al.l (|200G| ) to simulate the interaction 
between supernova ejecta and its companion. Generally, the 
radius of a companion decreases with metallicity, which is 
directly originated from the dependence of the periods of 
systems on metallicity (see last subsection). 

We see in these figures that Tycho G is consis- 
tent with our high- metallicity resul ts. The SG companion 
model used in [Marietta et al.l l|2000l) also matches with our 
high-metallicity resu l ts sin ce this model is obtained from 
iLi fc van den Heuvell (|l997l ). who used a method and a stel- 
lar evolution code similar to that in Paper I with Z = 0.02. 
However, the MS model in [Marietta et"al] (|200G| ) is obvi- 
ously departed from the (i?f^, ^/_Rf ^) plane for all metal- 
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Figure 8. Similar to Fig[5]but for Z = 0.0001. Top: acE = 1-0; 
Bottom: acE = 3.0. 



licities, which imphes that a solar model is not a reasonable 
model when simulating the interaction between supernova 
ejecta and its companion. 

There are also two islands in the (i?f^, j4/i?f^) plane 
with Z — 0.0001 and acE = 3.0 (see the bottom panel in 
FigH]). The reason is same to that interpreted in subsections 
[3T1 

We also showed the recurrent nova, U Sco, in these fig- 
ures. Since the orbital period of the system will reduce, the 
radius of the MS companion in the system should decrease, 
and_A/7?2^ increase based on the equation of lEggletonl 
l|l983t ). The position of U Sco will move towards right-lower 
region in (i?2^, j4/i?f^) plane, and may enter into the most 
probable region (see the bottom panel in Fig. [SJ • The bot- 
tom panel in Fig. [5] also shows that the recurrent nova, U 
Sco, probably explodes as a SN la. 



Figure 9. The distribution of the masses and the space ve- 
locities of companions in SNe la remnant for Z = 0.03 and 
OCE = 1-0. The cross represents the position of Tycho G, which 
is the potential candidate of the companion of Tycho's super- 
nova l lRuiz-Lapuente et aI.ll2004 lB ranc hl l2004) . and the length of 
the bars of the cross represents observational error. The position 
of the MS star of a recurrent nova, U Sco, is indicated by the 
filled square, where its space velocit y is its orbital velocity rela- 
tive to the mass center of the system llSchaefer fc Ringwaldlll995l: 
iHachisu et al...2000a.b.) . Top: Z = 0.03; Bottom: Z = 0.02. 



3.3 the masses and space velocities of companions 
in SNe la remnant 

As mentioned in section [1] in the single-degenerate 
model, supernova ejecta collides into the envelope 
of its companion after SN la explosion and strips 

some hydrogen-ri ch ma te rial fr om the surf ace of 

Wheeler et ai] Il975l: 



Frvxell & Arnett' 'l98l': 'Taam & Frvxell.' 


198d 




Liven. Tuchman & WhcUer 1992: 


200(\ 


Marietta et al.l (j2000|) ran several 



'1984'; 'Chugai 
Marietta ct al 



two-dimensional numerical simulation of the collision 
between supernova ejecta and its companion, where the 
companion is a MS star, a subgiant (SG) star, or a red giant 
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Figure 10. Similar to Fig[9]but for Z = 0.01 (Top) and Z = 0.004 
(Bottom). 

(RG) star. They found that the hydrogen-rich material of at 
least O.IS Mq is stripped from the en velope of a companion. 
However. [Meng. Chen fc Hani |20o3) used a simple analytic 
method while a mor e phy sical companion model than 
that of [Marietta et ahl (|2000l ) to simulate the collision and 
found that the stripped materi al may be a s low as 0.035 
Mb . Observationally, Mattila et al.l (|2005l ) and iLeonardI 
(12003) showed that the amount of stripped hydrogen-rich 
material is less than 0.02 Mq which can be neglected in 
comparison with companion mass. After the collision, the 
companion gains a kick velocity, whic h is much low er than 
the orbital velocity of the companion (iMariett a et al . 2000'; 
iMeng. Chen fc Han|[2007l ). 

For the reasons above, we assume that the mass of 
a companion is not changed by the collision of supernova 
ejecta and its space velocity after the collision is equal to 
its orbital velocity at the moment of explosion. Figs. [9lto ll2l 
present the distributions of the masses and the space veloci- 
ties of companions in SNe la remnant for different metallici- 
ties. Observationally, mass can be deduced from the spectral 
type of a star combining with the star's surface gravity and 
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Figure 11. Similar to Fig |9] but for Z = 0.001 (Top) and Z = 
0.0003 (Bottom). 

luminosity, and space velocity can be obtained from proper 
motion and radial velocity combining with distance. Then, 
our results can be compared with observations. The cross in 
the figures represents the position of Tycho G, which is a 
potential candidate of the companion of Tycho's supernova 
(iRuiz-Lapuente et al] |2004 iBranchI 12004 iHernandez et all 
I2OO9I ). It is interesting that the position of Tycho G is well 
consistent with our high-metallicity results, although Tycho 
G wa s suspected not to be the companio n of Tycho's super- 
nova (|Fuhrmannll2005l : llhara et al.ll2007l ). 

We see in the figures that the masses and the space 
velocity both increase with metallicity on average, which is 
derived from the dependence of the final secondary mass and 
the final period on metallicity as shown Figs. [T]to|4j 

The bottom panel in Fig. [T2l shows the {M^^ , V) plane 
for Z — 0.0001 and acE = 3.0. There are also two is- 
lands in the figure. The low-mass high- velocity one is from 
EAGB channel and the high-mass low- velocity island is from 
TPAGB channel. The reason is same to that interpreted in 
subsections 13.11 

In these figures, we also show the position of the MS 
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Figure 12. Similar to Fig[9]but for Z = 0.0001. Top: acE = 1-0; 
Bottom: acE = 3.0. 



star of the recurrent nova U Sco, where its space velocity is 
its orbital velocity relative to the mass center of the system. 
The fate of U Sco as a SN la is also clearly shown in the 
bottom panel of Fig. [9] 



Figure 13. The distribution of the gravities of companions for 
different metallicities with acE = 3.0, where the gravities are in 
cm/s^. 



o 

it 




Figure 14. The distribution of the opening angles of the holes 
in SNe la remnant for different metallicities with cicE = 3.0. 



3.4 The distribution of the gravities of 
companions 

The surface gravity of the companion in a SN la remnant 
is another parameter which can be directly obta ined from 
spectral observations (jRuiz-Lapuente et al.l [2004h . Surface 
gravity will be changed by the impact of supernova ejecta. 
As depicted in subsection 13.21 companions after the impact 
of supernova ejecta are processing into thermal equilibrium. 
Meanwhile, the thermal equilibrium of companions is still 
not reestablished at the moment of supernova explosion for 
pre-explosion mass transfer. Considering the insignificant 
change of companion mass, the surface gravity of a com- 
panion at the moment of explosion may represent a real one 
after the interaction to some extent, and then we could make 
a comparison between the distribution of the calculated sur- 



face gravities of companions and observations. We show the 
distributions of the surface gravities at the moment of ex- 
plosion for various Z and acE = 3.0 in Fig. 1131 uqe does not 
significantly change the distributions expect for Z — 0.0001. 
Double peaks also appear for the case of Z = 0.0001 in the 
figure. The low- gravity peak is from the TPAGB channel and 
the high-gravity one is from the EAGB chamiel. We see from 
the figure that the surface gravity for high- met allicity stars 
are generally larger than that for low-metallicity ones, which 
is directly derived from the larger masses and smaller radii 
of the companions with a high Z. Most stars have surface 
gravity \ogg in the range of 2-5. The potential companion 
of Tycho's supernov a, Tycho G, has a surface g ravity log g 
between 3.0 and 4.0 (|Ruiz-Lapuente et al.|[2004 ). Seen from 
the figure, it is difficult to judge the metallicity of Tycho G 
based on its surface gravity. 
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3.5 Percentage of SNe la with polarization 
spectrum 

After the interaction between supernova ejecta and its com- 
panion, the companion star carves out a conical hole of open- 
ing angle 30° - 40° in the supernova ejecta and the hole 
will never disappear because t he ejecta are moving super- 
sonically (iMarietta et al.|[2000l ). The aspheric configuration 
of supernova ejecta may r eveal itself by pol arization spec- 
trum near maximum light teasen et al ] i2004 ). With the ad- 
vance of spectropolarimetric observations, the nature of SN 
la asphericity becomes an importa nt relevant test of the sin - 
gle degenera t e prog enitor scenario ()Wang fc Wheeleill2008h . 
iKasen et al] (|2004l ') studied the effect of the hole on the 
spectra of SN la near maximum light and found that if the 
opening angle is larger than 20° , polarization spectrum may 
be obtained from any viewing angle. Considering the limit 
of spectropolarim etric observations (p olarization should be 
larger than 0.2%. [Leonard et alll2005l ). a polarization spec- 
trum can be detected when viewing angle is smaller than 
120° (defining that the viewing angle down the hole is 0°). 

The opening angle of the ho le, 6, in a SN la rem nant 
is mainly determined by A/R^^ l|Marietta et al.ll2000l '). We 
use 



e = 66.82 - 9.285(A/i?| 



- 0.3784(A/i?f'')' 



(6) 



to calculate the opening angle of the hole. The equation 
is fitted from the half-mass angle of th e distribution of 
stripped material given by , Marietta et aP (|2000,) and the er- 
ror of the equation is less than 1%. Note that the half-mass 
angle may be different from the opening angle and equa- 
tion (6) may overestimate the opening angle by about 5°. 
Here, we assume that equation (6) is valid for all metal- 
licities since metallicity may not significantly affect the 
interaction between supe rnova ejecta and its companion 
jMeng, Chen fc HanllioOTi ). 

Fig. [ij] show the distribution of the opening angles of 
the holes in SNe la remnants for different metallicities with 
acE = 3.0. The results of acE = 1.0 are similar to those 
of acE = 3.0 expect for the case of Z = 0.0001. Double 
peaks also appear in Fig. 1141 where the high-^ peak is from 
the TPAGB channel and the low-6 peak is from the EAGB 
channel. We see in the figure that the opening angles of 
the holes in all SNe la remnants are always larger than 20° 
even if the overestimate by equation (6) is considered. We 
assume that only when the view angle is smaller than 120° , 
polarization can be detected. Then, the percentage of SNe 
la with polarization spectrum is 1 — i~'^°°60 _ q 75 



4 DISCUSSIONS 
4.1 Tycho G 

It is believed that SN la is from the thermonuclear run- 
away of a CO WD in a binary system. The CO WD ac- 
cretes material from its companion to increase its mass. 
When its mass reaches its maximum stable mass, it 
explodes as a thermonuclear runaway and almost half 
of the WD m ass is converted into radioactive nickel-56 
l|Branchll200i ). Two progenitor models of SNe la have com- 
peted for about three decades. On e is a single-degenerat e 
model, which is widely accepted (jWhelan fc IbenI Il973h . 



In this model, a CO WD increases its mass by accret- 
ing hydorgen- or helium-rich matter from its companion, 
and explodes when its mass approaches the Chandrasekhar 
mass limit. The companion may be a m ain- sequence star 
(WD -f -MS) or a red-giant star (WD-I-RG') llYungelson et al] 



I995I: iLi fc van den Heu^ Il997l: 



Nomoto et al.l 



two channels 
ied and 



19991 : 
WD 



Langer et al. 



+ MS 
observations also 



Hachisu et alJ Il999al lb[ 
2OOOI ). Between the 
model is widely stud- 
the channel 



IS 

led and some obse rv ations also uphold tne cnanne 
jHachisu et alj Il999allbt iNomoto et all 1 19991: iLanger et al 



l2000l:lHan fc Podsiadlowskill2004l : Paper P. [Hachisu fc Kate 
(|2003al lb^ suggested that supersoft X-ray sources (SSSs) 
may be good candidates for the progenitors of SNe la, 
where some of SSSs be l ong to WD-I- MS channel. Recently , 
iHachisu fc Katol l|2005l . l2006al lbl) and iHachisu et all (|2007j ) 
showed that several recurrent novae is possibly the progen- 
itor of SNe la and some of them be l ong t o the WD-fRG 
channel. Observationally, iPatat et al.l |2007) suggested that 
the companion of the progenitor of SN 2006 X were an early 
RGB star. However, IHachisu et"ai] l|2008l ) argued a WD 
+ MS nature for this SN la. Considering a smaller calcu- 
lated Galactic birth rate of SNe la f rom WD -I- MS chan- 
nel th an that derived observationally (|Han fc Podsiadlowskil 
I2OO4I : Paper I) and some WD 4 - RG systems as the ca n- 
didates of SNe la progenitors (IHachisu fc Katol l2006al fbl: 
IHachisu et "ai]|2007l : IParthasarathv et aLll2007h . WD + RG 
channel should be carefully investigated further although 
some BPS results showed a small cont ribution of WD -|- 
RG channel to the birth rate of SN e la l|Yungelson fc Livid 
1 19981 : iHan fc Podsiadlowskil 120041 '). The other progenitor 
model of the SNe la i s a doubl e dege nerate model (DD, 
llben fc Tutuk^ 1 1984 IWebbinkI 1 19841 ). in which a sys- 
tem consisting of two CO WDs loses orbital angular mo- 
mentum by gravitational wave radiation and merges. The 
merger may explode if the total mass of the system ex- 
ceeds the Chandrasekhar ma ss li mit (see the revie ws by 
iHiUebrandt fc Niemeverl I2OO0I and lLeibundgutll2000l ). The 
birth rate calculated fr om this channel was comparable with 
the observational rate ('Han"l998'; 'Yungels on fc Lividll998l . 
2000; Tutukov fc Yungelson 2002). SN 20 03fg and SN 2005hj 
is likely the c a ses from the DD channel (jHowell et al.ll20()6l : 
iBranchI [200^ : lOuimbv. Hoflich fc Wheeled |2007| ). In addi- 
tion, KPD 1930-1-2752 may be an excellent candidate of DD 
SN la progenitor, whose total mass (~ 1.52Mq) exceeds 
the Chandrasekhar mass limit and whose orbital shrink- 
age caused by gravitational wave radiation will lead to the 
merger of the bi nary in about 200 Myr, much smaller than 
the Hubble ti me (ICeier et al"]|2007l). However , pleas e pay at- 
tention that lErg ma. Fedorova fc YungelsonI l|200l[ ) argued 
that, from detailed binary evolution calculation, the final 
mass of KPD 1930-1-2752 is smaller than the Chandrasekhar 
mass limit due to a large amount of mass loss during evolu- 
tion. 

A good way of discriminating between the many SN la 
progenitor scenarios is to search the companion of a SN la 
in its remnant. Unless the companion is another WD (DD 
channel, in which it has been destroyed by the mass-transfer 
process itself before explosion), it survives and shows some 
special properties in its spectra, whic h is originated from 
the contamination of super n ova ejecta ([Marietta et alj|2000l : 
iRuiz-Lapuente et al.l |2004| : iBranchI l2004h . Tvcho's super- 
nova, which is one of only two SNe la observed in our Galaxy, 
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provides an opportunity to address obs ervationally the iden- 
tificat ion of the surviving companion. iRuiz-Lapuente et al.l 
searched the region of the remnant of Tycho's super- 
nova and suggested that Tycho G, a sun-l ike star, is the com - 
panion of Tycho's supernova. Although, llhara et al.l (|2007| ) 
argued that the spectrum of Tycho G does not show any 
special properties, which seems to exclude the possibility of 
Tycho G to be the companion of Tycho' supernova, the anal- 
ysis of the chemical abundances of the Tycho G upholds th e 
companion nature of the Tycho G l|Hernandez et alll2009l) . 

Interestingly, some integral properties of Tycho G, i.e. 
the mass, the space velocity, the radius and the surface grav- 
ity, are all consistent with our binary population synthe- 
sis results (see Figs. [5] [9] and I13p . Then, Tycho G is very 
likely to be the companio n of Tycho's supernova (see also 
iMeng. Yang fc Gengj 120091 ') . If Tycho G were the compan- 
ion of Tycho's supernova as shown by our BPS results, it 
would challenge one's understanding about the physics of 
SNe la, such as the interaction between SNe la ejecta and 
companions. 



4.2 The simulation of the interaction between 
SNe la ejecta and companions 

In section [1] we have shown that [Marietta et al.l ()2000h ran 
several high-resolution two-dimensional numerical simula- 
tions of the collision between the ejecta and the compan- 
ion. They claimed that about 0.15 — O.17A/0 of hydrogen- 
rich material is stripped from a MS or a SG companion. 
They also found that in a sense of the collision, there is 
no difference whatever the companion is a MS star or a 
SG star. However, th eir results did not o btain confirma - 
tion by observation s |Mattila et al.l l2005l : iLeonar dl l2007l ). 
iMeng. Chen fc Hani l|2007^ used a simple analytic method 
to simulate the the interaction between SNe la ejecta and 
companions and found that for a given condition, more 
hydrogen-rich material is stripped from the envelope of a 
SG companion than that of a MS companion. To discuss 
the validity of their method, they repeated their work us- 
ing the s ame method in their paper and the companion 
model in [Marietta etHI l|2000l ). and found that the result 
of SG model is similar to that in [Marietta et"al] l|2000D . 
while the amount of the stripped hydroge n-rich material 
of MS model is much smaller than that in [Marietta et al.l 
ll2000t). The reason of the diffe rence is that the SG model i n 
[Marietta et al.l ([2000l ) is from [Li fc van den Heuvell l|l997l ). 
who used a method and a stellar evolution code similar to 
that i n Paper I. However, the MS model in Marietta et al.[ 
l|2000l ) i s a solar model, which is no t a typical case (see Figs 
[5lto[5|). lMeng. Chen fc Hani (|2007[ ') claimed that the differ- 
ence of the results bet ween their analytic m ethod and the 
numerical simulation in [Marietta et all l|2000l ) is mainly de- 
rived from the different stellar structure of the companion. 
Thermal equilibrium is not reestablished for the compan- 
ion star at the moment of supernova explosion since the 
mass transfe r is still processing be fore supernova explosion. 
Addit ionally, [Marietta et al.1 1|200 0^ and ' Meng. Chen fc Han! 
l|2007l ) argued that the luminosity of a companion after the 
impact of SNe la ejecta would rise sharply to about 5000 Lq, 
which is too high to compare with that of Tycho G. There- 
fore, considering that the properties of Tycho G is consistent 
with our BPS results, a detailed numerical simulation about 



the interaction between supernova ejecta and its companion 
should be en couraged by a more p hysical companion model 
than that in [Marietta ct al.' (2000). 

Recently, .Pakmor c t al. (20oA used a more physical 
co mpanion model and s imilar numerical simulation to that 
in [Marietta et al.l (|2000l ) to recalculate the interaction be- 
tween supernova eje cta and companion. T hey found a simi- 
lar results to that in [Marietta et af] (|2000l ). and at the same 
time they claimed that under some special situation, re- 
sults consistent with observation may be obtained, and then 
they claimed that theor y does not conf li ct wit h observation. 
Based on the resuh in [Pakmor etHI (|2008[ '), if A/R2^ is 
smaller than 6, the amount of striped mass from compan- 
ions should be larger than 0.02 A'Iq. According to the results 
in this paper (see figures 5 to 8), A/R2^ is always smaller 
than 6, and then the range of the stripped mass is from 
0.07 Mp) to 0.1 6 M(7), whic h is consistent with the res ults of 
[Marietta et al.[ ([200oi) and [Meng. Chen fc Had (|2007[ ). The 
effort of [Pakmor et all (|2008l ) then do not overcome the 
confl i ction between theory and observations (jMattila et al.l 
[20051 : II eonar d| [2007l '). The impact of a SN la on its compan- 
ions still should be studied carefull y. Maybe, unsymm etrical 
explosion plays an important role (jPlewa et al.[[2003 ). 



4.3 Polarization of SNe la 

As an important diagnostic tool for discriminating among 
SN la progenitor systems and theories of the explosion 
physics, spectropolarimetry provides the direct probe of 
early-time SN geometry. The essential idea is as follows: 
electron scattering dominates a hot young SN atmosphere 
and its nature is highly polarizing. For an unresolved 
source with a spherical distribution of scattering electrons, 
the directional components of the electric vectors of the 
scattered photons cancel exactly, yielding zero net lin- 
ear polarization. An incomplete cancellation will be de- 
rived from any asymmetry in the distribution of the scat- 
tering electrons, or of absorbing material overlying the 
electron-sca ttering atmosphere. Then, a net p olarization 
is res ulted ()Leonard fc Fihppenkd [20051 : [Wang fc Wheeierl 
[2008h . Single-degenerate model provides a natural way to 
produce the asymmetry. The exist of a companion may 
change the configuration of supernova ejecta and a polariza- 
ti on spectrum is e xpected. In this paper, w e use the results 
of [Marietta et al.[ 0OOa ) and [Kasen et al.[ ( (2004.) to calcu- 
late the percentage of SNe la with polarization spectrum 
and found that about 75% of all SN la may be detected by 
spectropolarimetry. However, this result critically depends 
on the following assumptions: (i) all SNe la are from single- 
degenerate progenitor systems, (ii) [Marietta et al.[ (|2000l ) 
showed reliable simulations in the sense that a hole is in- 
deed for med and does not quickly close with time and (iii) 
Kas en et al. (2004) provided a reasonable simulation of po- 
larization spectrum resulted from the exist of a hole in a 
supernova ejecta. 

It is likely that the single-degenerate model is only one 
of the reliable models, such as the prompt componen t in the 
two-co mp o nent model sugge s ted by iScannapieco fc BildstenI 
l|2005h and [Mannucci et al.l (|2006f ). At present, any defini- 
tive conclusion about DD model is premature, and this sce- 
nario can naturally result in an asymmetry of distribution 
of supernova ejecta. One mechanism is the rapid rotation 
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of a WD before supernova explosion, which leads to the 
change of the stellar shape. Another one is that there may 
be a thick accretion disk around CO WD and the disk is an 
origin of the asymmetry of the configuration of supern ova 
e.iecta (see the rev iews bv lHillebrandt fc Niemevej|2000l and 
lLeibundgutll200ol for details about DD model). Additionally, 
explosion mechanism itself may also prod uce the asymmetry 
and a pola rization spectrum is expected (jPlewa et al.ll2004l : 
iKcisen &r Plewa 2005i) . Then, it is probable that the percent- 
age of SNe la with polarization spectrum estimated in this 
paper is a lower limit. At present, almost all SNe la, which 
are observed by sp ectropolarimetry, ha d various degrees of 
polarization signal (jLeonard et alll2005l ). 



4.4 The percentage of SN 1991T-like 

SN 1991T is an overluminous event and has a rather broad 
light curve (Amis = 0.95 ± 0.05, where Amis is the mag- 
nitude difference b etween its maximum and 15 days later 
IPhillips et aLllTgggl ) , which is often taken as an indication of 
a large '^^Ni mass llHoflichl et al.1 1 19951 : iNugent et all 1 19971 : 
iPinto fc Eastmanll200ll ). Its spectrum also showed some spe- 
cial properties, i.e. dominated by Fell and Felll lines at 
maximum light, while the spectrum of a normal SN la is 
dominated by Sill line. A large ^^Ni mass may well ex- 
plain the peculiar s pectral appear ance (Ijefferv et aL 1 19921 : 
iMazzaH et"aLl 119951 ). Nevertheless. iKasen et al.1 (|2004 ) sug- 
gested a second, physically very different route to explain 
the spectral peculiarities of SN 1991T — one could be peer- 
ing down an ejecta hole. They found that their synthesis 
spectrum down the ejecta hole can be comparable with that 
of SN 1991T. We calculate the frequency of the special event 
for different metallicities by assuming that if view angle is 
smaller than the opening angle of the hole in a SN la ejecta, 
the SN la show properties of SN 1991T-like. The opening 
angel is from FigfH The calculated birth rate of SN 1991T- 
like supernovae is from 13% to 14% for different metallici- 
ties. The rate slightly i ncrease with Z and is a little larger 
than th e estim atio n by iKasen et al.1 (|2004l ') (~ 12%). Both 
iBranchI (|200ll ) and iLi et al. 1 l|200ll ) gave the observed rate 
of SN 1991T-like from 3% to 5%. If taking SN 1999aa-hke 
as SN 1991T-like events , the birth rat e of SN 1991T/SN 
1999aa-like i s 20% ± 7% llLi et al.ll200ll '). Our resuhs match 
with that of iLi et al] (|2001 ) within errors. So, it is possi- 
ble that SNe 1991T-like have not any special properties in 
physics except for the viewing angle of an observer. 



5 SUMMARY AND CONCLUSION 



Incorporating the results of iMeng. Chen fc Hani ||2009D into 
Hurley's rapid binary evolution code, we have carried out 
a series of binary population synthesis calculation and sys- 
temically study the properties of the companions of SNe la 
for different Z at the moment of explosions. We give the 
distributions of the masses, Mi^, the radii, i?!^, of com- 
panions and the periods, Psn and the ratios of separations 
to radii, A/R2^, of final binary systems for various Z at 
the moment of supernova explosion and find that generally, 
increases and -Rf^ andPsN decrease with Z, while the 
distributions of A/R2^ are similar for all metallicities. These 



parameters can be applied to constrain the numerical sim- 
ulation of the interaction between the ejecta of a supernova 
and its companion. In addition, these parameters can help 
to judge whether a WD + MS system may explode as a SN 
la or not. We also show the distributions of some integral 
properties, i.e. the masses, the radii, the surface gravities 
and the space velocities of companions for different Z after 
the interaction. The distributions can provide help to search 
companion in supernova remnant. Especially, some integral 
properties of Tycho G (a potential candidate of the com- 
panion of Tycho's supernova), such as mass, radius, surface 
gravity and space velocity, well match with our BPS re- 
sults. This fact may challenge our understanding about the 
physics of SNe la, especially the interaction between super- 
no va ejecta and it s comp anion . Using the r e sults simulated 
by [Marietta et al.l l|2000l ') and iKasen et al.1 (120041 ') . we find 
that about 75% of all supernovae can be detected by spec- 
tropolarimetric observations. If considering that there may 
be different progenitor models of SNe la and different explo- 
sion mechanisms, the percentage could increase. 
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